EXPONENTIAL MARTINGALES AND TIME INTEGRALS OF 

BROWNIAN MOTION 



VICTOR GOODMAN AND KYOUNGHEE KIM 



Abstract. We find a simple expression for the probability density of 
f exp(ZJ s — s/2)ds in terms of its distribution function and the distribution 
function for the time integral of exp(B s + s/2). The relation is obtained with 
a change of measure argument where expectations over events determined by 
the time integral are replaced by expectations over the entire probability space. 
We develop precise information concerning the lower tail probabilities for these 
random variables as well as for time integrals of geometric Brownian motion 
with arbitrary constant drift. In particular, E[ exp [6/ f exp(B s )ds)] is finite 
iff 8 < 2. We present a new formula for the price of an Asian call option. 
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1. Introduction 

Time integrals of one-dimensional geometric Brownian motion have appeared in 
financial models where certain expected values are the computed prices of 'Asian 
options'. Approximate values for some of these expectations were obtained in RSJ. 
Some fundamental work in |Y(92)| and |GY| focuses on the distribution and density 
functions of time integrals. Matsumoto and Yor [MY 1-3] studied time-integrals of 
the form 



/ exp(— 2B S — fis)ds 
Jo 



for [i > and obtained some interesting identities in law (Theorem 1.1 and 2.3 
in |MY3| ) . In a recent paper Bertoin and Yor [BY) considered the time- integral 
of Brownian functions over an infinite time interval. The distributions of time 
integrals over an infinite interval are known [D2j : 



,oo 1 

/ exp(— 2B S — fj,s)ds — 

Jo 2 7 M; 



where the law of 7 M / 2 is a gamma distribution with index /x/2. 

We present new relationships between the density functions and distribution 
functions for these random variables. Our results provide a precise description 
of the lower tail of these distributions and we settle several moment questions 
involving an exponent which is the reciprocal of a time integral. We also consider 
other moments where the exponent includes further exponential terms involving 
Brownian motion. Certain of these double exponential terms have a finite expected 
value. 

l 
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We let Bt denote a one-dimensional Brownian motion. Our notation for a time 
integral of exponential Brownian motion is somewhat nonstandard. We let Mt 
denote the simple exponential martingale 



M t = exp(B t - -) 



and we define its time integral as 



At = I MM. 



o 

It is known (see Dufresne, Corollary 2.3, [D]) that for a > and (i > \ 

^ exp ' < oo 
a + pA t 

We show that this moment, in fact, is finite for (3 = |. 

In [Y92,02] and BTW] the authors use A t to denote the time integral of expo- 
nential Brownian motion without drift, or an integral with a drift other than —1/2. 
We obtain some distributional properties of the random variable At . In particular 
we present a formula for its probability density and we obtain a sharp lower tail 
estimate for its distribution. 

Our method of argument uses a delicate change of measure, where the process 

(1-1) S t :=B t -21og(l-|^ t ) 

becomes a standard Brownian motion. Here, y is a positive constant. Of course, 
this path translation has a singularity at the random time defined by 

(1-2) A T ^ = - 

y 

Therefore, we will consider sample paths only for t strictly less than t^. Matsumoto 
and Yor (Section 2 in [MY3] ) considered a change of measure where 

=B t W -21og(l + 7^ t ) 

and the random coefficient 7 M is independent of the Brownian motion and has a 
gamma distribution with index \x. In contrast, our process B t is defined without 
this random factor and our process exists only up to the explosion time Too. 



2. The Change of Measure 



We follow the Girsanov formalism: A sufficient condition for a process 

B t :=B t - [ 6 s ds 



to be a Brownian motion over a compact time interval [0, T] is that the process 



(2.1) A t :=exp( / 8 s dB - | / 9 2 s ds) 



t x rt 
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be a martingale over the time interval (see [KS]). Bt is a standard Brownian motion 
w.r.t. the measure Q defined by 

(2.2) § = A T 

In view of equation (| 1 . 1 [) , we consider as our choice of 9t the process 

(2.3) *_ 2 <^_ 

Notice that Rq = —y and that R t is defined up to the random time Too. The process 
R t has the following convenient and remarkable property: 

Lemma 2.1. The process given in equation h2.S\) satisfies the SDE 

(2.4) dR t = R t dB t - \R\dt 
for t<T oa . 

Proof. This is a simple calculation using the fact that dMt — MtdBt. □ 



Remark. By writing equation \2.1$ in its integral form 



t 1 rt 



R t = -y+ / R s dB - - / R z s ds 



o 2 j 



we see that R t is essentially the exponent of the Girsanov density process it gener- 
ates. This unusual property of Rt allows us to analyze the behavior of A t through 
a change of measure. 



Definition 2.2. For each n = 1, 2, . . . let r„ denote the stopping time given by 

T n = va£{t : R t < — n} 

Although each stopping time, and Too as well, depends on the choice of y , we will 
omit mentioning their dependence on this parameter unless we explicitly change 
the value of y. We see from equation (|1.2|) and the fact that Rt is negative that 

T n < Too and also lim r n = t^ a. s. 

n — >oo 

Proposition 2.3. For each n = 1, 2, . . . the process 
(2-5) A { t n) = exp[ y + i? tAr „ ] 

forms a martingale. Moreover, the process 

B f -21og(l-|A AT J 
is a standard Brownian motion for t < T w.r.t. the probability measure 

dQ = k^dP 
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Proof. We choose a Girsanov density process as in equation (|2.1|) by setting 

6 S = i? s l{ s<Tii } 

With this choice, 8 S is a bounded adapted process and hence satisfies a Novikov 
condition (see Corollary 5.13 of JKS]). The Novikov condition is sufficent for A t to 
be a martingale. In our case, 

At = exp ( / R s dB - - / R 2 s ds) 
Jo z Jo 

It follows from Lemma |2. II that the exponent above is precisely Rt/\r n + U- That is, 

A t = exp(y + iW n ) 

This proves the martingale assertion of the proposition. Moreover, the calcula- 
tion in equation (|2.3p shows that 

/ Ms = 21og(l-|^ tAT „) 

and so the Girsanov theorem implies that the process 

B f -21og(l-|v4 Mr ,J 

is a standard Brownian motion on compact time intervals with the change of mea- 
sure given by At- 

□ 



3. The Correspondence Between B t and B t 

Definition 3.1. For fixed n — 1,2, .. . and y > we let B t denote the process 
appearing in Provosition \2.3[ That is 

(3.1) S t := J B t -21og(l-|A tAr J 
We also define 

M t = exp(S t - t/2) 

and 

A t = f M s ds 
Jo 

We note that all quantities in this definition depend on our choice for n and y. 
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Proposition 3.2. 

(3 - 2) M - = (1 -iw 



(3-3) 1 + V -A tATn - ' 



2 IAT " l-|A tAr „ 
and 



(3-4) iW„ = 



y- 1 + \A tKTn 



Proof. From Definition 13.11 we have 



M t = exp(B t -t/2-2 log(l - |^ tAT J) 



(l-|A tA rJ 2 

Now if £ < t„ then 

-(1 + -A) = 2 ' - 2(1 ipp = 



Hence 



V 7 1 

i + -A t = 



2 1 1-lAt 
Finally, for t in this same range, 

M t 



Rt = - 



y- 1 - hA t 



2 

M t (l - \A t f 

y- 1 \At 

= -yM t {l- V -A t ) 

yM 



1 + lAt 



2 ' 

These equalities hold up to the time t„ , and these are the assertions (|3.3|) and (|3.4j) 
in the proposition. □ 
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Proposition 3.3. If f(x, z) is a nonnegative Borel-measurable function and y > 
then 



(3.5) 



E[f{M u R t ) ■ A t <2/y] 



E[f( 



-M, 



■)exp(- 



Proof. For fixed n and y > we consider 

E[f(M t ,R t ) ; Tn >t] 

Now 

f(M t , R t )l{ Tn >t} = f(M u Rt) exp{-R t - y) exp{R t + y)l {Tn>t} 

Since this function vanishes for t > r„, each time parameter may be replaced by 
t A r„. This allows us to apply Proposition 12.31 where we take 

A t n)l {r„>t} = exp(i? t + 2/)l{ T „> t } 
We obtain the identity 

(3.6) E[f(M t ,Rt) ; r„ X] = B Q [/(Af t , i? t ) exp(-i? t - 2/)l {r „ >i} ] 



Each term in the r.h. expected value can be expressed in terms of the Brownian 
motion B t . We use the identities in Proposition 13.21 to see that 

/(M t ,i? t )exp(-i? t -y)l {Tri>t} = 

/(M t (l - \A t )\- M * )exp( M * - y)l {Tn>t} 
1 V 1 + %A t y 1 + ±A t 

(l + §A t ) 2 + y^ + lA 

Moreover, the event r„ > t equals the event 

mini? s > — n 

s<t 

which in turn equals the event 

Ms 



— < n 



The r.h. expected value in equation (|3.6[) is then 

... ^ M t . . M t . , 
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But since the integrand is nonnegative we may take the limit as n — + oo and 
obtain the limiting value 



Q[/( (TT|Xj2'~^ r TpI ) exp( ^Tp: - v)] 

In addition, since 



lim t„ = 

n — >oo 

as Too is defined in equation (|1.2|) we see that the limit of the l.h. side of equation 
(HOI) is 



E[f{M u R t ) ■ A t <2/y] 
This establishes the identity (13. 5|) of the proposition. 

□ 



Remark. Proposition \3.3\ is quite similar to Theorem 1 of [WH . The authors 
study Girsanov density processes and develop necessary and sufficient conditions for 
a Girsanov process to be a martingale. Theorem 1 shows that, in great generality, 
the expected value of a Girsanov density equals the tail probability for a certain 
stopping time. 

Our arguments proving our proposition are similar to those in [WHJ. In our 
case, we extend their result to include expected values of a function of the Brownian 
motion process multiplied by a Girsanov density. The choice f{x, z) = 1 is a special 
case of the theorem in [WH] . The reader may see this by making the choice 

W a + A t 

as required by Theorem 1. To define the correct stopping time, one should consider 
the process 

4 1 
Y{u)du := -21og(l - -A t ) 
a 

One may verify directly that Y(t) satisfies the functional equation mentioned in 
Proposition 1 of |WHj . 



4. The Distribution of A* 



Theorem 4.1. For a > the distribution of At is given by 



(4.1) Pr{ ^ t < a} = e -| i? [ exp (J^)] 

a + At 







Moreover, the random variable A t has a continuous, positive probability density 

Al 

M t 



function gt{a) which is simply related to the distribution functions of A t and jjf-: 
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ft 



(4-2) fft (a) = -Pr{A t <a}--Pr{-f < a} 



Proof. The first assertion of the Theorem follows from Proposition 13.31 by making 
the simple choice 

f(x,z) = 1 

Identity (|3.5[) becomes in this case 

Pr{A t < 2/y} = ff[exp( M \ - y)] 

= £[exp( ™\ )]e~y 

y +^A t 

This identity has the rather surprising corollary that the expected value above is 
finite. The integrand involves a double exponential of Browian motion. Using a 
result of Dufresne, Corollary 2.3 in [D], one can easily show that for (3 > h 

but here we have shown finiteness for (i = | . 

Notice that the integrand is a monotone function of y. If y varies over some 
positive interval 

(2/0,2/0 

then each integrand is dominated by the integrable random variable 

exp( ^f r rp; ) 

If a sequence {yk} converges to y in this interval, we apply the dominated conver- 
gence theorem to prove that the expected value converges to its value for y. Hence, 
the expected value is a continuous function of y. It follows that the distribution 
function of A t is continuous and we may write the identity as 

(4.3) Vr{A t < 2/y} = ff[exp( M \ )] e ~« 

V + 2 A t 

This establishes equation (|4.ip . 

We show the probability density exists by proving that the right hand expression 
in equation (|4.1|1 is differentiable w.r.t. a. To see this we consider the case of 
Proposition 13.31 for 

f(x, z)=x 

Identity (|3.5|) becomes 

Mt Mt 
E[ M,;A, < m . B|p - | ^ p exp(- TT V i -)]e-« 
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As in the previous case, the expected value 



is necessarily finite and again the integrand is monotone in y. Therefore, the ex- 
pression is a continuous function of y as we argued in proving (|4.1[) . However, this 
same expected value arises by formally differentiating the expected value 

(4-4) £[exp( M \ )] 

y 1 + %A t 

on the r.h. side of (|4.ip w.r.t. y. Now, as y varies over some positive interval 
(yoi Di) each integrand 

-2 M t M t , 

y /-i , i , ^ ex p( — i , i , ) 



it 

is dominated by the integrable random variable 

- 2 M t , M t 

y ° ^r 1 + 1^) 2 exp V 1 + 1 V 

So, the y— integral from yo to y\ 1 which equals 

, Mt , , Mt v 

is dominated by the product of y\ — yo and an integrable function. Therefore, the 
difference quotient for the expected value in (|4.4p converges as yo — > y\ and the 
limit is the entire expression for 

(4.5) -E[M t ;A t <2/y]e y 

The same argument applies to the case y± — ► yo so that expression (|4.4[) has a 
derivative which equals the expression (|4.5[) . Since the distribution function is the 
product of e~ v and expression (14. 4p we conclude that the distribution function of 
At has a continuous probability density. We differentiate terms in the identity ()4.3j) 
to obtain 

(4.6) - 2y- 2 g t (2/y) = - Pv{A t < 2/y} + E[M t ;A t < 2/y] 
Now, the expected value 

E[M t ;A t < 2/y] 

is another distribution function. Using the change of measure induced by the factor 
Mt, we see that 

B s - s = W t 

is a standard Brownian motion and so the expected value equals 



Pr{ / exp(W s +s-s/2)ds<2/y} 
Jo 
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We substitute this expression into equation ()4.6|) to obtain 

g t (a) = ^(Pr{A t < a} - Pr{ f cxp(B s + s/2)ds < a}) 
a Jo 

We see that the expression for gt(a) is strictly positive since 

f exp(B s + s/2)ds 
Jo 

is strictly larger than the random variable A t for each sample path. Finally we note 
that the random variable has the same distribution as the time integral above. 



4r= f e MBs -B t + t/2 - s/2)ds 

M t Jo 



has the same distribution as 



f exp(W t _ s + 1/2 - s/2)ds 
Jo 



where W s denotes a standard Brownian motion; we may change variables in the 
time integral to obtain the time integral of geometric Brownian motion with positive 
drift. This establishes assertion (|4.2[) of the theorem. 

□ 



Remark. We can rewrite the density formula {4 .Sty by combining the two probabil- 
ities. The density equals 



(4.7) 



g t (a) = \Pr{[ exp(B s - s/2)ds < a < [ cxp(B s + s/2)ds} 
a Jo Jo 



This expresses the probability density for A t in terms of a single condition on Brow- 
nian motion sample paths up to time t. 

The existence of a continuous probability density for A t is a corollary of Propo- 
sition 2 of Yor [Y(92)| . In the proposition a conditional density for A t is given as 
an integral transform of various transcendental functions. No explicit connection 
is made between the density and the distributions of At and At /Mt . 

In Dufresne [D] nice formulas are obtained for the density of a reciprocal of a 
time integral for some values of a drift parameter in the Brownian motion. Our 
choice corresponds to the choice of \i = —1, and density is given as an integral 
transform in [DJ. 



Remark. The density <?t(a) is a solution of the PDE derived in [BTW, equation 
(26)] where it is shown that time integrals of more general geometric Brownian 
motions have smooth densities. In particular, the random variable 

A_ 

has a smooth density since it has the same distribution as the time integral of 
exp(_B s + s/2). The PDE, equation (26), has a simple form in the case of the 
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density of At: 

dg d 2 , a 2 g , dg 

di = da 2 ^^' ~ da 
But, equation \4-*fy shows that a 2 g/2 is the difference of two distribution functions. 
Therefore, the PDE becomes 

dg d , d A t , dg 
at da 1 da M t da 

d 2 r A t 
oa 2 M t 

That is, the time derivative of the density is obtained by differentiating the density 
ofAt/Mt. 



Corollary 4.2. For each y > the process 

is a supermartingale. In particular, 

(4.8) E\Z t \ = e y Pr{A t < 2/y} 

so that E[Z t ] is a strictly decreasing function oft. 



Proof. The process 

Yt - M > 



y- 1 + \A t 

also satisfies the SDE that appears in Lemma [2~T1 A simple calculation shows that 

(4.9) dY t = Y t dBt - ~Y t 2 dt 

Therefore, the remark following Lemma I!Q1 applies to the process Y t : 
Since 

Y t - y = f Y s dB -If Y 2 ds, 
Jo 1 Jo 

Y t — y is the exponent of the Girsanov density process which Y t generates. A simple 
stopping time argument, similar to one in the proof of Proposition 12.31 shows that 

z t = cx P (r t ) 

is a positive local martingale. However, the integral identity (|4.ip of the Theorem 
implies that 

£[exp(Y t )] 

is decaying function of t. Consequently, Z t is a local martingale but it is not a 
martingale. □ 
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Remark. The process Yt appears implicitly in Lemma 2.1 of |BTWj . In order 
to derive PDE's for certain expected values involving At, the authors compute a 
diffusion equation for processes slightly more general than 

(Y t -\M t ). 

The PDE identities do not apply to equation |^.ff| ) since Z t is not a homogeneous 
function ofYf 1 . 

One can derive the corollary from Theorem 1 of [WHj . but the argument here 
connects the result directly to the behavior of time integrals of geometric Brownian 
motion. 



Corollary 4.3. 

(4.10) £[exp(i-)] = oo 

Proof. For each a > 

/ 2 s 

E[exp (-J-) ;A t <a}> exp(2/a) Pr{A t < a} 
At 

and equation (|4.1|) of the theorem implies that the r.h. expression equals 



E[ex ^^TAt )] 

But, this quantity increases as a — > and therefore the random variable 

is not integrable. □ 



5. Geometric Brownian Motion with Drift 



Definition 5.1. For each i/elue let A^ denote the time integral 

(5.1) A\ v) = / exp(B s + us - s/2)ds 

Jo 

The random variable A t in the previous sections is A^ with this notation. 
Theorem 5.2. For a > the distribution of A^ is given by 

(5.2) Pr{A C> < a } = fe-iEKa + A^ exp ( 2 ^ + "~ 

a + A\> 
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Proof. For y > we apply Proposition 13.31 to evaluate 

E[(M t r ; At < 2/y] 
The choice of f(x, z) = x" in the proposition gives the expected value 

t?i M t , M t 

£ W 2 ^- exp( ^TFr y)] 

Next, we multiply these expected values by exp(j/i/2 — v 2 t/2) so that 

M"exp(vt/2 - vH/2) = exp(i/B t - v 2 t/2) 

We use this exponential martingale factor to change measure in each integral so 
that the process 

B s = B s — vs 
is a standard Brownian motion for s < t. We see that 

E[M" exjp(vt/2 - v 2 t/2) ; A t < 2/y] = Pt{A[ v) < 2/y} 
while the other expected value equals 



E[(l + y -4^)- 2u eX p( eXp(i3t + " t 7; /2) 
This establishes identity (|5.2[) . 



2/)] 



□ 



Corollary 5.3. As a J, i/ie function 

exp(2/a) ft/a) ^ fl} 



increases. 



Proof. For the case that ^ = 1/2, the formula in (|5.2[) for the distribution function 
becomes 

(5.3) Pr^ 1 / 2 ) < «} = oe"* 25[ exp )] 

The integrand of the expected value in (|5.3|) increases as a [ 0. 

□ 



Corollary 5.4. The following expected value is infinite. 
(5-4) £[exp(-^— -)]=oo 
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Proof. The corollary states that -E[exp ( , 2 /2) )] = oo. Let F(a) denote the dist 



bution function of A^ 2 ^ and consider 



^xp(-^);4 1/2) <i] 

A t 

We use a standard argument that justifies integration by parts: 

e« = / -^re^ax 
Ja x 2 

so that 

(\ pi pOO rj 

£[ ex p (-jrm) i 4 1/2) ^ !] = / / 3« f « f («) 

' ' J0 Ja x 

f°° 2 2 /' a:A1 
= / — e* / dF(a)dx 
Jo x J Q 

oo o 

Z 2 

—e*F(x A ljax 



But, Corollary 15.31 implies that 

1 * , . fc 

—cFyx) > — 

a; 

for x < 1 where the constant A: > 0. Hence, the integral is infinite. 



6. Finite Exponential Moments 



Lemma 6.1. For any t > 

(6-1) E[ex P ( t * )]< 

2 J Q exp(B s )ds 



Proof. The expected value in (J6JJ) is 

^t exp (^p))] 

with the notation of Section 5. And, an upper bound for the expected value is 



f°° 1 
e + ^ Pr{exp ( ^(i/a) ) ^ ^ 

Let a; = e s to obtain the following form of the integral above: 
(6-2) 



p oo 1 

: / p r { — > 4 1/2) }e s ds 
Ji 1 2s " ( J 
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It suffices to prove that the integral from k to oo is finite where k is chosen so that 
on the interval of integration 

4 

- < * 

s 

For any s in the interval we have 

A[ 1/2) = I exp(B u )du > / exp(B u )du 
Jo Jo 

f 4/s 4 s f 4/s 

/ exp(B u )du > - exp(— / B u du) 
Jo s 4 J 



4 s 
s ' 4 



by Jensen's inequality. Since the random variable 

r 4/s 



B u du 



/3s 



is normal with standard deviation 



we may replace it with 



3s 

where Z denotes a standard normal random variable. The integrand in equation 
(|6.2p is dominated by 

14 / 2 \ 
Pr{— > -exp -=Z }e s 



2s ~ s 



'2>s 



< Pr{-2 > -^=Z}e s 
V3s 

since log(l/8) < —2. We may write this as 

Pr{-V3s > Z}e s 
Hence the integral in equation (|6.2[) . 



is dominated by 



/OO -i 
Pr{ - > 4 V2) }e s ds, 

e k + / Pr{-V3s" > Z}e s ds 

Jk 

k r , 3s n SJ 

i e + / exp( — — )e as < cxd 

Jk 2 



□ 
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Theorem 6.2. For any 9 < 2 and to > the process 

,0M t 

is a supermartingale for t > to • In particular, 

(6.3) £[ exp (— i)] <oo 

At 

and the expected value is a strictly decreasing function oft. 

Proof. We have seen (proof of Theorem l4.ip that A t /M t has the same distribution 

as 

I exp(B s + s/2)ds 
Jo 

Since this random variable is larger than A^ 2 \ Lemma 16-11 implies that 

£[exp(^)]<oo 

Now let 

e ct Mt 

(6.4) U, - 

for any fixed c > 0. If to > is sufficiently small, so that e ct °/4 < 1/2, we have 

£[exp(t/ to )] < oo 

We define t\ by 

e ctl /4 = (9 

and we claim that 

E[exp(U tl )] < oo 

The integrabliity of exp (^4 ) will follow because the choice of c is arbitrary. We 
first consider the SDE for Ut- 

e ct M t 

dU t = cU t dt + U t dB t - - — j^Mtdt 
4 Aj 



UtdBt + cU t dt - -^Ufdt 



4 

-Ctl 



= UtdB t + U t {c - Ae~ ct U t }dt 
Next, we compute the SDE for the process exp (Ut). 

dexp(U t ) = U t exp(U t )dB t + exp(U t )U t {c - 4e~ ct U t }dt + i exp(U t )U?dt 
In order to compute an expected value, we introduce the stopping times 

t„ = inf{i >t :U t > n} 
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and we obtain 



E[exp(U TnMl )] - E[exp(U tQ )] 

(■TnAtl 



= E[ / exp(U t )U t {c + -U t - 4e- ct U t } dt] 

J to z 

< E[ / exp(Ut)Ut{c+-U t - Ae~ ctl U t }dt\ 

/•T„Atl 1 

(6.5) = E[ / exp(U t )U t {c + ~U t - B^Ut} dt] 



Notice that if 
where 



'to 

U t >b 
b:=c{9^-\)-' 



then the integrand in equation (|6.5p is negative. It follows that 

^[exp(i7 r „ Atl )]-^[exp(i7 t0 )] 

<E[[ n 1 eMUt)U t {c+\u t - e^Utjl^^dt] 

Now we take the limit as n — > oo and, using the condition that Ut < b in the 
integrand, we see that 

E[exp(U tl )] < oo 



This establishes that the expected value in equation (|6.3|) is finite. It is a strictly 



decreasing function of t because the random variable -A- has the same distribution 



as 

4 1} which is a strictly increasing function of t. It remains to establish the 
supermartingale property. Corollary 14.21 implies that for any a < 1 the process 

\ exp ( — \ 

is a supermartingale because it is a concave function of a supermartingale. Now as 
a — > the pointwise limit of this process is 

|exp(— )) =exp(— — ) 

That is, the process in the statement of the theorem is integrable and is the limit 
of non-negative supermartingales. Hence, it is also a supermartingale. □ 



Remark. It follows from Theorem 1 5. j| that the distribution function of has the 
form 

a 2 exp(—2/a)K a 



where the function K a increases as a J. 0. So, if 9 > 2 in equation \6.3\) . the 
expected value is infinite. It is unclear for the case 9 = 2 if the expected value is 
finite. 
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Corollary 6.3. For any 9 < 2 and t > 



(6.6) E[exp(- t — — )]< 

J exp(_B s )ds 



Proof. Since the expected value in (|6.6|) is a decreasing function of t, it suffices to 

Mt 



show the expected value is finite for t arbitrarily small. Since ^ has the same 



distribution as 

f exp(B s + s/2)ds, 
Jo 

Theorem 16.21 implies that for any 9 < 2 

E[exp ( — )] < 

J exp(i3;5 + s/2)ds 

For a given 9 < 2, choose t so small that 



'exp(^) <2 



so that 



Mexp ( -. )] < oo 

V e-*/ 2 / exp(Sa + «/2)d« 



We see that the expected value above is larger than 

E[e*p( t ~ — - )] 
J exp{B s )ds 

and so this expected value is finite. □ 



Remark. Some related work on exponential moments of A t is mentioned in Yor 
[Y(92)| . Equation (l.e) of the article states that 

E[ 7T t exp ( " Ya) ] = ^TWt exp ( " ^ sinh " u ^ 

However, a difference in notation requires that A t be expressed with our notation 
as jA^ 1 ^ . By writing t' = t/4 we see that the formula in Equation (l.e) is an 
expression for 

E l i exp ( j— rj 



Equation \6. 6)) implies that the expected value is finite for all complex values of u 
such that Re(u 2 ) > — 1 and is analytic on this region. Since 



sinh 1 u = log (it + yl + u 2 ) 

one sees that the right hand expression also has an analytic extension. The formula 
has a singularity at the value u — i which corresponds to the infinite exponential 
moment of Corollary \5.4\ 



GEOMETRIC BROWNIAN MOTION 
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In addition, Theorem 4-1 of [D] provides an integral formula for the density of 
A^ 2 ^ and one can show that the expected value is finite (for 9 < 2) using Theorem 
4-1- The corollary provides a probabilistic proof of this fact. The result that the 
moment is infinite for 9 — 2 is new. 



An expected value considered in Y(92)| , |RSj , and in other works concerning 



time integrals in financial mathematics is the 'price' for an Asian call option. We 
present a new expected value for the simplest option price and indicate how to 
derive a corresponding formula for arbitrary constant drift and volatility values. 

Proposition 6.4. For any a > 

9 M 9 

(6.7) E[(A t -a)+ }=t-a + a 2 E[(a + A t )~ l exp )] 

a + A t a' 



Proof. The notation (A t — a) + involves the indicator function of the event A t > a 
so it is natural to consider 

E[ A t - 2/y ; A t < 2/y } 
= - 2 -E[\- V -A t ■ A t < 2/y] 

y 2 

We apply Proposition 13.31 where we make the choice 

f(x,z) = — 

Then 



-2f(M t ,R t ) = -2Mt 



M, 



y 2 

Equation (|3.5p shows that the expected value over A t < 2/y is given by 
— E[(l + f A t - 1 exp « - y)} 

On the other hand, 

2 2 
E[A t --]=t-- 

y y 

so we take the difference of these expected values and let | = a to obtain the 
result. □ 



Remark. The proposition requires the Brownian motion to have drift —1/2. To 
derive an expectation formula for arbitrary drift one can apply Proposition \3.3\ to 
the expected value 

E[{M t Y( A t - 2/y) ; A t < 2/y } 
The factor of (Mt) v is relevant for a change of measure (see the proof of Theorem 
\5.2\) so that the time integral will contain an arbitrary drift. To incorporate a 



20 



VICTOR GOODMAN AND KYOUNGHEE KIM 



volatility factor in the Brownian motion, one can make a simple time scale change. 
We let s = cr 2 s' so that 

-i f exp(B s )ds = { exp(B a 2 s i)ds' 
° Jo Jo 

This random variable has the same distribution as the time integral o/exp(cr_B s ). 



References 

[BY] J. Bcrtoin and M. Yor: Exponential Junctionals of Levy processes, Probabilty Surveys 
2, 191-212 (2005). 

[BTW] R. Bhattacharya, E.Thomann, and E. Waymire: A note on the distribution of integrals 
of geometric Brownian motion, Stat, and Prob. Letters 55, 187-192 (2001). 

[D] D. Dufresne: The integral of geometric Brownian motion, Adv. in Appl. Probab. 33, 

223-241 (2001). 

[D2] D. Dufresne: The distribution of a perpetuity, with applications to risk theory and pen- 
sion funding, Scand. Actuar. I. 1-2, 39-79 (1990). 

[GY] H. Geman and M. Yor: Asian Options, Bessel Processes and Perpetuities, Math. Finance 
2, 349-375 (1993). 

[KS] I. Karatzas and S.Shreve: Brownian Motion and Stochastic Calculus, Springer- Verlag 
New York (1991). 

[K] K. Kim: Moment Generating function of the inverse of integral of geometric Brownian 

Motion, Proc. Amer. Math. Soc. 132, 2753-2759 (2004) 
[MY1] H. Matsumoto and M. Yor: An analogue of Pitman's 2M - X theorem for exponential 

Wiener junctionals, Part I: A time-inversion approach., Nagoya Math. J. 159, 125-166 

(2000). 

[MY2] H. Matsumoto and M. Yor: An analogue of Pitman's 2M — X theorem for exponential 

Wiener functionals, Part II: The role of the generalized inverse Gaussian laws., Nagoya 

Math. J. 162, 65-86 (2001). 
[MY3] H. Matsumoto and M. Yor: A relationship between Brownian motions with opposite 

drifts via certain enlargments of the Brownian filtration, Osaka Journal of Mathematics 

38, 383-398 (2001). 

[RS] L.C.G. Rogers and Z. Shi: The value of an Asian option, J. Appl Appl. Probab. 32, 
1077-1088 (1995). 

[WH] B. Wong and C.C. Heyde: On the martingale property of stochastic exponentials, J. 

Appl. Probab. 41, 654-664 (2004). 
[Y(92)] M. Yor: On some exponential functionals of Brownian motion, Adv. in Appl. Probab. 

24, 509-531 (1992). 

[Y01] M. Yor: Exponential Functionals of Brownian Motion and Related Processes., Springcr- 
Verlag, New York. (2001). 

Mathematics Department, Indiana University, Bloomington, IN 47405 
E-mail address: goodmanvSindiana.edu 

Mathematics Department, Florida State University, Tallahassee, FL 32306-4510 
E-mail address: kimamath.fsu.edu 



